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ABSTRACT 

Using the observed rate of short-duration gamma-ray bursts (GRBs) it is possible to make predic- 
tions for the detectable rate of compact binary coalescences in gravitational-wave detectors. These 
estimates rely crucially on the growing consensus that short gamma-ray bursts are associated with 
the merger of two neutron stars or a neutron star and a black hole, but otherwise make no assump- 
tions beyond the observed rate of short GRBs. In particular, our results do not assume coincident 
gravitational wave and electromagnetic observations. We show that the non-detection of mergers in 
the existing LIGO/ Virgo data constrains the progenitor masses and beaming angles of gamma-ray 
bursts (e.g., 6j > 4° for M to tai > 20 M©, for uniform component mass), although these limits are fully 
consistent with existing expectations. We make predictions for the rate of events in future networks of 
gravitational- wave observatories, finding that the first detection of a NS-NS binary coalescence associ- 
ated with the progenitors of short GRBs is likely to happen within the first 16 months of observation, 
even in the case of a modest network of observatories (e.g., only LIGO-Hanford and LIGO-Livingston) 
operating at modest sensitivities (e.g., advanced LIGO design sensitivity, but without signal recycling 
mirrors), and assuming a conservative distribution of beaming angles (e.g. all GRBs beamed with 
9j = 30°). Less conservative assumptions reduce the waiting time until first detection to a period of 
weeks to months. Alternatively, the compact binary coalescence model of short GRBs can be ruled out 
if a binary is not seen within the first two years of operation of a LIGO-Hanford, LIGO-Livingston, 
and Virgo network at advanced design sensitivity. We also demonstrate that the rate of GRB triggered 
sources is less than the rate of untriggered events if 6j < 30°, independent of the noise curve, network 
configuration, and observed GRB rate. Thus the first detection in GWs of a binary GRB progenitor 
is unlikely to be associated with the observation of a GRB. 



1. INTRODUCTION 

The LIGO and Virgo collaborations have recently re- 
leased results from roughly a half year of observations, 
investigating the gravitational w ave (GW) sky at un - 
precedented levels of sensitivity (lAbadie et al.l l2012bh . 
They did not identify any gravitational wave sources, and 
thereby established new upper limits on the rates of a va- 
riety of p ossible GW events in the nearby (< 200 Mpc) 
Universe (|Abadie et al.|[201l|) . One of the most promis- 
ing sources for GWs detectable by these ground-based 
observatories is the coalescence and merger of a compact 
binary system: two neutron stars (NS), two black holes 
(BH), or one of each. 

There has been an active program of observing 
gamma-ray bursts (GRBs), focusing on rapid follow- 
up to determine afterglows and identify host galax- 
ies ( Soder berg et al . 2006; Panaitescu 2006; B erger et al.l 
120071: IPerlev et al.l \200m . As a result, there is grow- 
ing evidence that short/hard gamma-ray bursts are as- 
sociated with the mergers of either t wo neutron stars , 
or a neutron star with a black h ole (jFong et al.l 120101 : 
IChurch et all [201lt IBergerl l20ll . This consensus is 
based on noting that the physical timescales are com- 
mensurate, the short GRBs do not appear to be asso- 
ciated with star formation (and therefore are unlikely 
to be associated with supernovae), and the GRBs oc- 
cur far from the centers of their host galaxies. These 
studies have also provided rcdshifts for a subsample of 
short GRBs, thereby providing preliminary estimates 



for the rate densities of these events ()Nakar et al.l 120061 : 
[DieT3[20Tl . There is great interest in gravitational 
wave/electro magnetic multi-messenger observations of 
these GRBs (IMetzger & BergerllMa lEvans etaH\201% 
iBriggs et al.ll2012l ). as such systems would help confirm 
the first detections of GWs, elucidate the properties of 
GRBs, and potentially provide interesting measurements 
of the H ubble constant and the dark energy equation- 
of-state (ISchutzl [19861: IHolz fc Hughei [200a [Dalai et al.l 
12001 iNissanke et al.H2oTol )~ 

One of the most important properties of GRBs is the 
beaming of the gamma rays. This beaming directly re- 
lates to the total electromagnetic energy of the explo- 
sion, as well as the intrinsic event rate of the sources 
(as compared to the observed rate, which is a function 
of the ones that happen to point at us). Recent ob- 
servations of a jet break in the short-duration gamma- 
ray burst GRB 1110 20A suggests a b eaming opening 
angle of 6 3 ~ 3-8° (|Fong et al.ll20Ll . Other GRBs 
(e.g., GRB 070714B, GRB 070724A, and GRB 071227) 
have been found with beaming angles in the range 1— 
30° (|Fong et al.l [2011 ICoward et aL|[2012T ). while non- 
detection of a jet break in the light curve of GRB 
050724A places a lower limit on th e beaming of that 
burst of 9 



les, 



3 > 25° dGrupe et al.ll2006D . N umerical stud- 
on the other hand, find 0, < 30° ( Popha m et alj 



1991 IRosswog fc Ramirez-Ruiz|[200l IJanka et all 120061; 
Rezzolla et al.H2011l) . 

In this paper we estimate the limits that arise on 
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the beaming of short-duration GRBs based on the 
non-detection of GWs from associated binary systems 
in the recent LIGO/Virgo science run. We also make 
projections for the detection rate of binary systems, 
as a function of mass and beaming angle, for future 
networks of GW observatories. We take a conserva- 
tive lower limit on the observed rate density of local 
short GRBs of K G rb = lOyr^Gp c" 3 dNakar et al.l 
[200l |Diet3 l20m iCoward et all 12012 ), based primarily 
on BATSE and Swift observations. We emphasize 
that this rate is determined purely through observa- 
tions, although it is broadly consiste nt with the rates 
arisin g from population sy n thesis ( Belczvnski et al 
20061: lO'Shaughnessv et all 120081: jBelczvnski et al 



2010t lAbadie et al.l 120101: iDominik et all 120121) . The 



gravitational-wave limits presented here are based on 
observed GRB rates, and are therefore independent of, 
and complementary to, estimates based on population 
synthesis modeling. 

We assume that all short GRBs are associated with 
low-mass compact binary coalescence. There is de- 
veloping evidence that this is the case, with per- 
haps a small sample of nearby GRBs occurring from 
other m echanisms, such as flares from soft g amma re- 
peaters (jLevan et alJl2008t lAbbott et al.ll2008[ ). While it 
is conceivable that not all short GRBs are the result of 
binary coalescences, it is perhaps even more likely that 
not all binary coalescences result in GRBs. We thus ex- 
pect that our limits on the minimum beaming angle in 
Sec. [2] are low, and our estimates of the maximum wait 
time in Sec. [3] are high. 

2. LIGO S6/ VIRGO VSR2 

From July 2009 to October 2010 the LIGO and Virgo 
observatories conducted a search (S6/VSR2-3) for co m- 
pact binary coalescences (jAbadie et al.l l201ll l2012aH bl). 
They did not detect any gravitational-wave events, and 
thereby established upper limits o n the event rates o f 
coalescences in the local Universe (jAbadie et al.ll201ll ). 
The LIGO instruments operating during S6 were the 
4 km laser interferometers at Hanford, WA [H] and Liv- 
ingston, LA [L], while the Virgo [V] results were from a 
single 3 km laser interferometer in Cascina, Italy. The 
S6/VSR2 runs consisted of 0.09 years of HLV coincident 
data, 0.17 years of HL, 0.10 years of HV, and 0.07 years 
of LVQ 

We have taken the representative sensitivities pre- 
sented in Fig. 1 of lAbadie et al.l (|2011l ), and calculated 
the corresponding horizon distances, Rq, for H, L, and 
V, where horizon distance is defined to be the distance 
at which a given signal-to-noise (SNR), p, is measured 
for an optimally oriented (face-on ) and optimally loc ated 
(directly overhead) binary. From (Dalai et al.l ((20061) : 

R = AA^T 7 /p, (1) 



where A = ^5/96tt- 2 / 3 (GM/c 3 ) 5 / 6 c, and where the 
binary chirp mass is given by A4 = (™im2) 3 / 5 /( m i + 
tt^) 1 /' 5 . In this paper we are only interested in nearby 
sources (z < 0.2), and for simplicity neglect the redshift 

1 Due to compromised s ensitivity from the installation of an 
incorrect mirror, we follow Abadic et al. (2011) in removing the 
Virgo SR3 data from our analysis. 
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Fig. 1. — Maximum rate density of binary coalescences as a func- 
tion of total mass of the binary, given the non-detection of any bi- 
nary systems in the LIGO/Virgo S6/V2 observing runs. The dot- 
ted blue curve is the result from Fig. 4 ofEbadic et al. ( 20fJJ), set- 
ting an upper limit on the rate density, where each mass bin is aver- 
aged over a uniform distribution of component masses. The dashed 
green curve shows our result for this curve, where we have assumed 
a fixed SNR network threshold of p = 9.4. The lower sensitivity 
of the Abadie curve at higher mass is due to non-stationary noise. 
We calculate an effective SNR network threshold, as a function of 
mass, to match the Abadie curve, and we use this to calculate the 
equivalent curve assume that mi = 1.4 Mq (i.e., that the NS mass 
is fixed, and the black hole mass is given by mi = M tota i — 1.4). 



dependence of the chirp mass (IHolz fc Huehesl 12001 . 
The characteristics of the detector are encapsulated in 



h = 



J hi. 



-7/3 



s h (f) 



df. 



(2) 



with Sh(f) the noise spectral density of the detector. 
We follow the ap proach of LIGO's com pact binary coa- 
lescence searches (jAbadie et al.ll20l2bf ). and take /i ow = 
40 Hz for the LIGO detectors, and f\ ow = 50 Hz for 
Virgo, with upper limits set by the frequency of the in- 
nermost stable circular orbit. To calculate the waveform 
we assume the members of the binary are non-spinning, 
and make use of the stationary pha se approximation 
(jDalal et al1[200l lAbadie et al.ll2012bf). 

We follow the approach of iSchutzl (|20~llT) to combine 
the antenna patterns of the three different interferome- 
ters (HL, HV, LV, and HLV), taking into account the 
differing horizon distances (which are a function of the 
masses of the binary) as well as power patterns. The 
networ k weigh t ed ant enna power pattern, from eqs. (15)- 
(20) of lSchutzl ((2011 . is given by: 



P 



N\?,<P) = 
N 

{(Fl k (8A^)+Fl k (9,<p^))p^ k Rl k ) 



(3) 



k=l 

where F® k ( 



,4>,ip) and F 2 k (0, </>, ip) are the antenna 
patterns for the fcth detector, and p m in,k an d i?o,fe are 
the SNR threshold and horizon distances of that de- 
tector, respectively. The detection distance for a net- 
work is related to the antenna pattern: R(0, <jy) = 

y/ Pn{Q, >t>) I Pmiai where p m ; n is the network signal-to- 
noise threshold. The sensitivity also depends upon the 
orientation of the binary; integrating over all orientations 
results in a factor of 0.29 decrease in the mean detectable 



GRB beaming and GW detection 



3 



volume (jSathvaprakash fc Schut3l2009t ISchutzll2011l) : 



V 



0.29 1 

3 Pmin 



(4) 



Since the LIGO/Virgo network did not detect any 
gravitational wave sources during its last science run, we 
can combine all network configurations and correspond- 
ing coincident observational times to estimate a 90% up- 
per limit to the rate density: TZ = 2.3/(^VK; x Ati), 
where the sum is over the different detector networks 
configurations, Ati is the amount of observational time 
for network configuration i, and the factor of 2.3 is in 
accordance with a Poisson process (see the discussion 
below eq. [6]). We plot our results in Fig. [1] assuming 
a combined network thresh old of p = 9 . 4. Th is is to 
be compared with Fig. 4 of lAbadie et al.1 (j 2 1 lh . which 
calculates the same quantity through detailed analysis 
of the GW data stream from the LIGO and Virgo in- 
struments ([Brady et al.ll2004r i. We have tuned our SNR 
threshold to agree for low mass binaries, but our re- 
sults begin to deviate at higher mass (M to tai ^ 10 Mq), 
since the signal shifts to lower frequencies and is there- 
fore more sensitive to non-stationary noise in the detec- 
tors (in part because there are fewer cycles to integrate 
against). To get a sense of the importance of this effect, 
we calculate the effective network SNR threshold which 
we would need to apply, as a function of mass, to match 
the rate l imits which come o ut of the full analysis pre- 
sented in lAbadie et al.1 (|2011[ ). From a value of p = 9.4 
at Mtotai = 3.5 Mq, the non-stationary noise degrades 
the sensitivity of the instruments at higher masses, lead- 
ing to p = 10.7 at 11-14 M Q , and p = 11.1 at 20- 
25 Mq. In what follows we incorporate this mass de- 
pen dence into our effect ive SNR thresholds. We note 
that lAbadie et all (|2011l ) assume a uniform distribution 
of component masses for their binaries. We also con- 
sider the case where the neutron star is restricted to have 
mi = 1.4 Mq, and the mass of the companion is given 
by m 2 = Mtotai — m i- Because this entails higher mass 
ratios for higher mass binaries, it decreases the overall 
gravitational-wave strength of the sources in comparison 
to the uniform distributions, and therefore decreases the 
detectable volume. As can be seen in Fig. [TJ this re- 
sults in a negligible effect at low mass, but rises to a 
factor of 3.6 in the rate at Mtotai = 15 M Q and 6.8 at 

Mtotai = 25 Mq. 

We define the beaming angle, 9j, to be the half open- 
ing angle of one of the two polar jets of a gamma-ray 
burst. The fraction of the sky, covered by the beamed 
gamma rays is given by /b = 1 — cos 9j . Given the 
paucity of data on the beaming of short GRBs, it is 
premature to assume knowledge of the distribution of 
beaming angles. We therefore will assume that all short 
GRBs have a fixed beaming angle, 6j, with the under- 
standing that this fixed angle provides the same results as 
the appropriate average of the true distribution of beam- 
ing angles. In other words, 1/(1 — cosOj) = JP(9)/(1 — 
cosd)d8, where P(9) is the true distribution of beam- 
ing angles. Assuming all short GRBs have compact bi- 
nary progenitors, the implied rate density of these coa- 
lescences is given by TZ = 7?. G rb / fb = TZ G ^ B / ( 1 - cos 0j ) . 
If only a fraction, /cbc , of short GRBs result from com- 
pact binary coalescence, then the rate density becomes 
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as a function of 

the total mass of the binaries. We plot results for three different 
rates of short GRBs, 72. = 3, 10, and 30yr Gpc , and for two 
different distributions (mi is uniform from 1 to iWtottJi and mi = 
1.4 Mq). Given the lack of detection of binaries in LIGO/Virgo 
S6/V2, the beaming angles of short GRBs will be greater than 
what is plot ted in 90% of cases. We use a network threshold which 
matches the Abadie ct al. (2011) results (see text for details). 

TZ = /cbc^grb/(1 - cos0j). In Fig. H we plot the 90% 
lower limit on the beaming angle as a function of the 
mass of the progenitors. We take the observed rate of 
short GRBs to be 3, 10, or SOyr^Gpc" 3 dNakar et~aTl 
l2lMlGuetta fc Stellal[200llCoward et alJl2012D . For ex- 
ample, a rate of 3yr _1 Gpc~ 3 can be thought of as a 
very conservative estimate of TZ — 2yr _1 Gpc~ 3 with a 
very conservative estimate of /cbc = 0.5. This sim- 
ple analysis suggests that models of short GRBs with 
progenitors of mass M to t a i > 20 Mq (uniformly dis- 
tributed in component mass) and with a beaming angle 
6j < 4° are inconsistent with existing LIGO /Virgo data. 
This weakens significantly for more realistic masses and 
mass ratios, with a minimum beaming angle of ~ 1° 
at Mtotai ~ 3Mq. These limits are completely consis- 
tent with observations and expectations. The current 
LIGO /Virgo data is on the verge of providing interesting 
astrophysical constraints, which suggests that the next 
generation of detectors should provide quick detections, 
or alternatively, the lack of quick detections would pro- 
vide strong lower limits on the beaming of short GRBs. 
We explore these constraints in the next section. 

3. ADVANCED LIGO/VIRGO 

We now calculate the expected detection rate of short 
GRB progenitors in the advanced LIGO and Virgo de- 
tectors, as well as additional detectors in Japan (KA- 
GRAfl and India (IndlGOfl. The advanced LIGO de- 
tectors are expected to begin operation in ~ 2015, and 
it is hoped that the LIGO and Virgo observatories will 
achieve their target advanced detector sensitivities by 
~ 2017, with the Japanese [J] and Indian [I] detectors 
operating at comparable sensitivities by ~ 2020. We 
assume an identical noise curve for each of these in- 
struments, given by the representative advanced LIGO 
noise curves in LIGO document TO9OO288-v30 with 

2 gweenter . icrr .u-tokyo . ac . jp/en/ 

3 www. gw-indigo . org/tiki-index . php 

4 dec . ligo . org/ cgi-bin/DocDB/ShowDocument?docid=2974 
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with the probability of waiting a time r before detect- 
ing the first event given by e _rA . We define i nr st as the 
waiting time by which, in 90% of cases, the first event 
will have been observed: ifirst = — ln(0.1)/A = 2.3/A. In 
Fig. [3] we plot tfirst as a function of the beaming angle, 
9j, for the HL and HLV networks. If one is interested in 
the waiting time by which the first event has been seen 
in 50% or 99% of the cases, the 90% waiting times are 
multiplied by 0.3 or 2, respectively. We have assumed 
/cbc x = 10yr _1 Gpc~ , and we have considered 

NS-NS equal mass binaries, with mi = mi — 1.4 Mq. 
We do not employ the mass-dependent threshold correc- 
tion factors from the S6/V2-3 analysis derived in the 
previous section, since it can be argued that the non- 
stationary noise is less likely to be a problem in the ad- 
vanced configuration of these detectors Q and, regardless, 
the functional form would not a priori be expected to 
match that of the lower sensitivity detectors. The rep- 
resentative curves in Fig. [3J can be rescaled to other pa- 
rameter values of interest. For example, if either T^grb 
or /cbc is down by a factor of 10, then the waiting 
times, ifi rs t, arc all multiplied by the same factor of 10. 
A change in the network configuration similarly results 
in an overall shift in the expected rates, and therefore 
an overall shift in ifirst' The relative waiting times for 
other networks, Tfirst, are presented in Table [TJ scaled 
to the HL (no SRM) value. The HLVJI network, with 
all detectors at the advanced zero-detuning high laser 
power sensitivity, has a waiting time that is a factor of 
0.15 that of the HL (no SRM) curve, independent of the 
beaming angle. If we consider larger mass systems, the 
waiting time is correspondingly shortened (see eqs. ([5]) 
and For example, if we consider NS-BH systems, 

with M NS = 1.4M and with M B h = 10 M and 20 M , 
then the waiting times until first detection are a factor 
of 0.11 and 0.05 shorter than those for the NS-NS case 
presented in Fig. [3J It is to be noted that our results are 
roughly consistent with the completely independent rate 
estimat es from population s ynthesis and observed binary 
pulsars (|Abadie et al.ll2010l) . 

Alternatively, the binary origin of short GRBs can be 
falsified (at the 90% Poisson confidence discussed above) 
if no coalescences are observed with a full network (HLV) 
operating at design sensitivity (zero-detuning high laser 
power) over a period of 3 years. This limit comes from 
taking 9j = 90°, which corresponds to none of the GRBs 
being beamed (which is already inconsistent with ob- 
servations). If we take a conservative upper limit of 
6j = 45°, we find that the binary origin can be falsified 
(at the 99% level) in 70 months for HL (no SRM), and 
in 22 months for HLV. However, even if short GRBs are 
not the result of binary mergers, we nonetheless expect 
a population of merging systems, and these should b e 
observable by future observatories (jAbadie et al.ll2010f ). 

The predicted event rates for the various networks are 
related to the 90% waiting times by a factor of 2.3. If an 
HL network, using conservative noise curves (no SRM), 
takes 16 months before a first detection assuming a bi- 
nary of mass IAMq-IAMq and a beaming angle of 
0j = 30° (see Table [T|) , then the predicted event rate 
for this network is 2.3/1.33 yr = 1.7yr -1 . These results 

6 Peter Shawhan, private communication. 



TABLE 1 

Mean detectable volume and wait times for the detection 
of binary coalescence associated with short grbs in 
future gw detector networks. the network snr is 
taken to be 10, and the volume is calculated for a 

1.4 M0-1. 4 Mq BINARY. Tfl rst IS THE WAITING TIME UNTIL FIRST 
DETECTION, SCALED TO THE VALUE FOR THE HL NETWORK WITH 
THE "NO SRM" NOISE CURVE. THE LAST THREE COLUMNS LIST 
THE 90% WAIT TIME FOR FIRST DETECTION (iN MONTHS) FOR 
THREE DIFFERENT VALUES OF THE BEAMING ANGLE. 



/low = 10 Hz. We take the target design sensitivity to be 
given by the ZER0_DET_high_P . txt curve, corresponding 
to zero-detuning of the signal recycling mirror, and high 
laser power. We also consider an early, less sensitive in- 
carnation of the detectors resulting from the absence of 
signal recycling mirrors, given by the N0_SRM.txt curve. 
It is possible that in 2015 the worldwide GW detector 
network will consist solely of HL in this lower sensitivity 
configuration. 

We calculate the mean detectable volume, V, of a vari- 
ety of ground-based networks. Our results are presented 
in Table [TJ where in all cases we have assumed a network 
SNR threshold of p = 10, and our sources are taken to be 
equal-mass binaries with mi = ni2 = 1.4 Mq. The mean 
detectable volume of a network is expected to scale as 
M 5/2 (see eqs. Q, (J3J), and (g])), although this relation 
is imperfect, since the scaling also depends on the shape 
of the noise curve. We fit the mass dependence to the 
functional form 

where Vo is the detection volume for a binary with 
M to t a i = 2.8 Mq, and we follow the previous section and 
fix mi = 1.4 M Q and m 2 = M tot ai - m 1 . The "no SRM" 
noise curve yields p = 1.39, while our fiducial advanced 
LIGO noise curve yields » = 1.30. These fits are good to 
20%atM total = 30M Q E| 

We are interested in how quickly networks of advanced 
ground-based gravitational-wave detectors can be ex- 
pected to see their first binary coalescences associated 
with short GRB progenitors. This rate is a function of 
the rate of observed GRBs (in gamma rays), 7£ G RB) the 
sensitivity and configuration of the detectors, the mass 
distribution of the GRB progenitors, and the beaming of 
the GRBs. The event rate of detectable binaries for a 
network of GW observatories is given by 

A = ^(M to tai)ft GRB /cBc/(l - cosOj). (6) 

How long will a given network have to wait before seeing 
its first event? This is described by a Poisson process, 

5 For completeness, we also mention results for the 
ZER0_DET_low_P.txt noise curve, corresponding to a lower laser 
power. We find V = 0.02 Gpc 3 for HL (no SRM), V = 0.034 Gpc 3 
for HLV, and p = 1.37. 
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Fig. 3. — Waiting time until first detection, ifl rs ti as a func- 
tion of the beaming angle, 9j. Results are shown for both un- 
triggcred and triggered GW observations of GRB progenitors, 
where we have assumed the observed local short GRB rate is 
^■GRB = 10yr~ 1 Gpc -3 . In 90% of cases the waiting time will 
be less than the values indicated. We plot results for two different 
GW networks: Hanford+Livingston, operating without a signal re- 
cycling mirror (a potential early version of advanced LIGO) and 
Hanford+Livingston+Virgo, operating at the design sensitivity (at 
high laser power). Even in the pessimistic case (HL, no SRM, all 
short GRBs are beamed with 8j = 30°), the first untriggered bi- 
nary detection is expected in less than 16 months. For more sub- 
stantial GW networks, the expected wait time may be less than a 
month (e.g., if 0j ~ 10°). The dotted lines show the wait times for 
GRB triggered GW observations, for a range of values of p. The 
GRB trigger provides a time and sky position, thereby reducing 
the required SNR threshold and increasing the detection rate be- 
cause the sources are assumed to be face-on. On the other hand, 
the triggering GRB rate is given by the 8j = 90° values, and is 
not enhanced by beaming. The dotted lines show the equivalence 
beaming angles, at which the rates of triggered and untriggered 
GW observations of GRBs match each other; this is found to hap- 
pen at 6j ~ 30° for p ~ 7.5. For lower values of 8j, untriggered 
observations occur more frequently than those triggered by GRBs, 
and the first observed binary GRB progenitor will be seen first 
(and perhaps only) in GWs. 

are consistent with those of iCoward et al.l (|2012f ). 

4. TRIGGERED VERSUS UNTRIGGERED 

An important question when considering the future 
GW detection of short GRBs is whether triggered or un- 
triggered detections are more likely. A short GRB trigger 
improves the sensitivity of the GW search by reducing 
the need to marginalize over all times and sky positions. 
In addition, because the gamma-rays are thought to be 
beamed, a GRB trigger is expected to be face-on, thereby 
increasing the signal in GWs over a source with a random 
inclination. In other words, the strongest GW emitters 
happen to also be the ones observable in gamma-rays. A 
GW network is therefore substantially more sensitive to 
a GRB triggered source. On the other hand, for small 
beaming angles the rate density of GRB progenitors in- 
creases (approaching oo as 9j —> 0). For sufficiently small 
9j the untriggered GW detection of a GRB progenitor 
will dominate over a triggered GRB, even considering 
the additional sensitivity in the latter case. We are in- 
terested in establishing whether triggered or untriggered 
GRBs arc more likely for upcoming networks of GW ob- 
servatories. 

We have calculated the waiting time and event rates 
for untriggered observations of GRB progenitor systems 
above. We now consider the equivalent calculation in 
the case of a GRB trigger. We assume that the GRB 
is face-on, which improves the sensitivity of the net- 



work by a factor of 1/0.29 (see eq. (jl])). We now esti- 
mate the reduction in SNR threshold due to the know n 
time and sky position of the source (|Dalal et al.1 120061) . 
We assume exp(— p 2 /2) oc l/# of templates, where in 
the untriggered case we took p = 10. If we take this 
threshold to have been based on roughly one year of 
observation, the existence of a GRB trigger now re- 
duces the observational window down to ~ 10 sec, for 
a reduction in the number of templates by a factor of 
~ 10 6 . If t he sky localization i n the untriggered case 
is ~ 5 deg (jNissanke et al.l 120 111) , then compared to a 
full-sky search (41, 253 deg ), the reduction in number of 
templates is a factor of ~ 10 3 . The total number of tem- 
plates is down by a factor of 10 9 , which for the equivalent 
false alarm rate would imply that the SNR threshold is 
reduced to p ~ - 21og(cxp(-10 2 /2) x 10 9 )) = 7.7. 

In Fig. |3] we compare the wait times for untriggered 
and triggered GRBs. We find that the wait time for a 
triggered GRB is equivalent to the untriggered case if the 
average GRB beaming angle is 9j ~ 30°. This result is 
independent of the network configuration, the individual 
noise curves, and the assumed assumed GRB rate, and is 
weakly dependent on the specific value of the threshold 
improvement due to the reduced number of templates. 
If the GRBs have an average beaming value of 9j = 20° , 
we would predict the rate of untriggered GRBs to be 
roughly double that of triggered ones. This increases 
to a factor of 10 if 6j = 10°. Alternatively, the rates 
are equal if 9j = 30°, and triggered GRBs are found at 
double the rate of untriggered ones if 9j = 45°. It is to 
be noted that this process can be inverted, and the wait 
time before first detection (and in between the first few 
detections) may be used to infer the beaming angle of 
GRBs. In addition, the relative rates of triggered and 
untriggered GRBs will help establish the beaming, and 
will be an important test of consistency when compared 
with explicit determinations of the beaming distribution 
based on GW measurements of the inclination of GRB 
sources. 

As discussed in Sec. [TJ recent observations have mea- 
sured short GRBs with 9j ~ 10°, indicating that it is 
likely that the rate of untriggered GRBs will be greater 
than the rate of triggered ones, and implying that the 
first detection of a binary system which is a progeni- 
tor of a short GRB will not be triggered by a GRB. 
Although triggered GRBs may be less frequent than un- 
triggered ones, multi- messenger observations of these sys- 
tems holds tremendous s cientific potential, and should 
be aggressively pursued (jBloom et al.l l2009h . Further- 
more, the increase in psychological confidence of detec- 
tion given coincident observation may play a large role 
in the initial detections. It is to be emphasized that 
our triggered rates assume the existence of an all-sky 
short GRB monitor operating contemporaneously with 
advanced GW networks. 

5. DISCUSSION 

We have explored the connection between the ob- 
served short GRB rate, the beaming angle of short 
GRBs, and the predicted rate of detectable binary sys- 
tems associated with progenitors of GRBs in networks of 
gravitational- wave observatories. 

We have shown that existing LIGO/ Virgo data pro- 
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vides preliminary constraints on the beaming angle and 
mass distribution of short GRB progenitor systems. For 
example, we find that short GRB progenitors of mass 
■Mtotal > 20 M Q (uniformly distributed in component 
mass) and with beaming angles of 6j < 4° are ruled out 
by existing LIGO/ Virgo data. These constraints, while 
novel, are fully consistent with our current understand- 
ing of the short GRB engine and rates. 

We have analyzed the observed rate of short GRB 
progenitors in future networks of GW detectors. We 
find that, even in the pessimistic case of only two de- 
tectors (HL) operating at conservative sensitivity (with- 
out a signal recycling mirror), in 90% of cases we would 
expect a first detection of a binary within 16 months 
if the GRBs are beamed within 6j = 30°, and within 
55 days if 6j = 10°. The expected event rates are 
1.7 yr" 1 (6j = 30°) and 15 yr" 1 {6, = 10°). We find that 
the HLV network, operating at zero-detuning and high 
laser power, would shorten these times to 4.9 months 
(9j = 30°) and 17 days (9j = 10°), with corresponding 
event rates of 5.6 yr -1 and 49 yr -1 . Alternatively, the 
binary coalescence model for short GRB progenitors can 
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be ruled out if an HLV network does not observe a binary 
within the first two years of observation. 

Finally, we have shown that the rate of GRB trig- 
gered observations of GW systems associated with GRBs 
is lower than the rate of untriggered observations if 
9j > 30°. This result is independent of network, noise 
curve, and GRB rate, and when coupled with recent ob- 
servations of small beaming angles for short GRBs, sug- 
gests that the first detections of GRB progenitors with 
advanced GW networks will not involve the observation 
of GRBs. 

We conclude that, assuming short GRBs are the result 
of the merger of compact objects, and assuming that the 
resulting gamma-rays are beamed, the first detection of 
gravitational-waves from binary coalescence associated 
with a GRB progenitor will be untriggered, and may oc- 
cur within months of operation of a modest network of 
ground-based gravitational wave observatories. 
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